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ORGANIZED TEAMWORK IN BASIC RESEARCH 


“What are the limitations of the organized teamwork method of research? Can it be 
effectively applied in carrying forward basic, fundamental, scientific inquiries? It has been 
questioned, for example, whether the teamwork method, which has unquestionably been 
proved effective and powerful in the type of applied research important to industry, can 
be expected to supersede the efforts of the'‘inspired, individual worker in the realm of pure 
scientific research. 

“The older generation of natural scientists—notably such British giants of the 19th century 
as Tyndall and Huxley—insisted that he who hopes to be successful in uncovering Nature's 
secrets must search in a humble spirit animated solely by a desire for truth and with no 
thought in mind of rewards or of the ultimate use, if any, to which the truths revealed may 
be put. . . . Clearly, the motive back of industrial research is quite different. . . . Yet this 
does not necessarily mean . . . that the methods which have proved so valuable in indus- 
trial research cannot be applied to research in pure science. . . . Motives differ, but methods 
may have values which are independent of motive. 

“On the one hand, the industrial researcher seeks answers to specific industrial problems 
and cannot ignore the fact that he is expected to show a reasonable probability of—must, 
indeed, in the long run, produce—financial returns commensurate with the costs of his 
searches. It is clear that few, except, perhaps, some of our major industries, can afford 
free-roving individual workers and gamble on the remote possibility that, sooner or later, 
they may turn up with a discovery which will be of value to that particular industry. 

“When it comes to writing a prescription for research, therefore, the doctors differ. There 
is the apparently incurable romantic group which insists that scientific truths come to us as 
special revelations granted to the very rare, inspired, self-sacrificing, and faithful scientific 
worker in a magical ‘flash of genius.’ There are others who hold to the statement, attributed 
to Edison, that genius is 1 per cent inspiration and 99 per cent perspiration. Important 
discoveries both in pure science and in engineering science and industry have been made 
by both methods. For example, Langmuir’s work with the General Electric Company in the 
development of the gas-filled electric lamp may be cited as an example of the value to 
industry of research which stemmed from a spirit of scientific curiosity, and began as an 
inquiry in pure science. Morgan’s epoch-making discoveries in genetics on the other hand, 
were made possible through organized teamwork, Industrial research, it is true, does aim 
to anticipate the possible discoveries of individuals through teamwork—but, here again, 
notable industrial advances—new products and processes—are still being developed by indi- 
vidual, independent inventors, although it seems certain that, as the equipment essential 
to productive research increases in scope and cost, the isolated inventor will face increasing 
difficulties. : 

“We are reminded, in this connection, of an observation by the British author Hilaire 
Belloc . . . ‘The conquests of physical science,’ he observes, ‘have been due to the minute 
and extensive observations conducted by vast numbers of men, and,’ he concludes, ‘there- 
fore, for the most part by the unintelligent’... . 

“In this statement Belloc was contrasting the remarkable result stemming from the work 
of many with the great-man theory of history—the idea embodied in the famous statement 
of Marshall Foch, that it was not an army that crossed the Alps, it was Hannibal. We may 
not like Belloc’s use of the word ‘unintelligent,’ but he put his finger on a truth which 
constitutes the major element of strength and power in science and engineering—the power 
of an ever-growing body of scientific and technical knowledge and understanding, and the 
fact that not only the individual genius but vast numbers of men of less than extraordinary gifts 
can add, through organized research, their bit to this ever-accumulating power. . . In science 
and engineering, thousands of workers, known and unknown, have piayed a part in building 
the huge accumulation of knowledge and experience from which the power of technology 
stems. It is precisely on this fact that modern organized and directed scientific and engineer- 
ing research relies.”—James Kir Fincu in Engineering and Western Civilisation, New York, McGraw- 
Hill Book Company, Inc., 1951, pp. 311-13. 


of Washington was established in December, 1917, to coordinate investi- 

gations in progress and to facilitate the development of engineering and 
industrial research in the University. Its purpose is to aid in the industrial 
development of the state and nation by scientific research and by furnishing 
information for the solution of engineering problems. 


The scope of the work is three-fold: 

(1) to investigate and publish information concerning engineering prob- 
lems of a more or less general nature that would be helpful in 
municipal, rural, and industrial affairs ; 

(2) to undertake extended research and to publish reports on engineer- 
ing and scientific problems ; 

(3) to provide opportunities for graduate engineers to conduct research 
under conditions that will most effectively prepare them for pro- 
fessional service. 


For administrative purposes the work of the Station is organized into nine 
divisions : 
. Aeronautical Engineering 
. Chemical Engineering, Industrial Chemistry 
. Civil Engineering 
. Electrical Engineering 
. Forest Products 
Geology 
. Mechanical Engineering 
. Mining, Metallurgical, Ceramic Engineering 
. Physics Standards and Tests 


The control of the Station is vested in a Station Board consisting of the 
President of the University, the Dean of the College of Engineering as 
chairman, the Director of the Station, and members of the faculty, representa- 
tive of the administrative divisions. The Board determines the character of 
the investigations to be undertaken and supervises the work. 

The Station offers a substantial number of research fellowships to highly 
qualified graduate students who work under the direct supervision of the 
faculty of the various divisions. Results of major investigations are published 
in the form of bulletins and reports. Current research findings from the 
Station’s fellowship projects, as well as accounts of its activities, appear in 
the quarterly journal, The Trend in Engineering at the University of Wash- 
ington. Reprints of articles by members of the engineering faculty and graduate 
students published in recognized technical journals are also issued by the 
Station. Requests for copies of the publications and inquiries for information 
on engineering and industrial problems should be addressed to the Director, 
Engineering Experiment Station, University of Washington, Seattle 5. 
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Fic. 1. X-RAY DIFFRACTION AND SPECTROSCOPIC UNIT 


B. S. R. Sastry, Research Assistant in the Engineering Experiment Station, is adjusting the diffractometer. In the center is the 
X-ray spectrograph, and at the right, the electronic control unit and recorder. 
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Analytical Aspects of X-Ray Diffraction 


James I. MUELLER 


Associate Professor of Ceramic Engineering 


The rapid advancement 
of the relatively new sci- 
ence of X rays has caused 
many of its potential ap- 
plications to be either over- 
looked or neglected. In the 
main, this is due to the fact 
that too few workers in the 
fields of pure and applied 
science are aware of the 
possibilities of this tool. 
Therefore a review of some 
of the fundamentals and a 
discussion of a few of the 
applications may be timely. 

The use of X-ray diffraction and fluorescent X-ray 
spectrographic techniques is becoming increasingly 
important in research and control in many branches 
of mineral industry. Applications in the fields of 
metallurgy, mining and prospecting, mineral prepara- 
tion, and ceramics have reached the point where, to- 
day, a modern, up-to-date diffraction laboratory is 
recognized as one of the necessary facilities for teach- 
ing in these fields. (See Fig. 1.) 


J. I. Mueller 


Fundamentals 


X rays are part of the electromagnetic spectrum 
and are similar to ordinary light rays in that they 
exhibit the properties of refraction, diffraction, and 
interference. A complete discourse on the generation 
and properties of X rays is not within the scope of 
this paper, which will be limited to a brief explana- 
tion of diffraction phenomena. 

Diffraction is actually part of the phenomenon of 
scattering. The formulation and explanation has been 
simplified by Bragg’s analogy with the reflection of 
ordinary light from a mirrored surface, i.e., consider- 
ing that the array of atoms within a crystal forms a 
series of parallel planes, the normal distance between 
these being noted by the symbol, d. (See Fig. 2.) 
Inasmuch as this d spacing depends upon the size of 
the atoms making up the crystals and the class in 
which the material crystallizes, it becomes evident 
that no two materials have an identical set of d 
spacings. This, together with the fact that certain 
extinction laws prohibit all planes from giving a 
reflection of the incident X-ray beam, makes it pos- 
sible to use the d spacing pattern as a “finger print” 
of the crystal. 
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The Bragg law is based on the fact that the path 
difference between the portion of the incident X-ray 
beam reflected by the successive planes must be equal 
to an integral number of wave lengths in order that 
the diffracted beam will be reinforced. (See Fig. 2.) 
This law is stated as follows: 


nd=2d sin 6 
where 
n=order of diffraction (1, 2, 3, etc.), 
\= wave length of the incident radiation, 
d=interplanar spacing, 
6=diffraction angle. 

In the application of X-ray spectroscopy still an- 
other phenomenon must be considered. When a pri- 
mary X-ray beam strikes an absorbing material with 
sufficient energy to eject an inner orbital electron 
from the atoms making up the absorber, a secondary 
X-ray photon, whose wave length is characteristic of 
the absorbing material, is emitted. By diffracting this 
secondary (fluorescent) X-ray beam by means of a 
single crystal, we can again apply the Bragg law, and 
by measuring the diffraction angle, 6, we can de- 
termine the wave lengths present in the spectrum. 
These can then be compared to known wave lengths 
to determine the composition of the material. 

The chief advantages of X-ray spectroscopy over 
optical spectroscopy are that (1) a much larger 
sample may be utilized, (2) the concentration of the 
material is unimportant, (3) the X-ray spectrum has 
fewer lines than the visible spectrum, and (4) 
analysis is independent of the physical state of the 
element. 

Uses of X-Ray Diffraction 
Qualitative Analysis 
Probably the greatest use of X-ray diffraction has 


been in the analysis of powder samples to determine 
the crystalline constituents present. The method is 
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based on the fact that each crystal has its own charac- 
teristic patterns of d spacings and, hence, of the posi- 
tion of the diffraction lines. In addition, the location 
and types of atoms within. the various planes cause 
the intensity of the diffracted lines to vary. It is then 
apparent that each crystalline material will have a 
set location of all its diffraction lines for a given wave 
length, and that each line will have a given relative 
intensity. The sample is usually finely ground (—200 
mesh) in order that the particles will have a high 
degree of random orientation. 

This randomness is necessary in order to compare 
patterns of the same material, since preferred orienta- 
tion would cause the relative intensity of certain lines 
to be lesser or greater than normal. Figure 3 shows 
a typical diffraction pattern of a sample containing 
two unknowns obtained in the Norelco Geiger- 
counter diffractometer in the X-Ray Diffraction Lab- 
oratory of the School of Mineral Engineering. 
Indexing of this pattern results in the data shown in 
Table I. A set of index cards has been published by 
the American Society for Testing Materials which 
gives the results of various diffractionists on more 
than 2,500 different crystalline materials. Each mate- 
rial is indexed according to the d value of the three 
most intense lines. Therefore, by selecting the strong- 
est lines of the pattern in Fig. 3 and checking the 
cards in the ASTM index, we find that the combined 
patterns from cards of sodium chloride and potas- 
sium chloride (Fig. 4) match the lines in Fig. 3. 

All samples to be studied may not be in powder 
form, e.g., metal specimens. The same procedure 
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Fic. 3. X-RAY DIFFRACTION PATTERN OF SODIUM CHLORIDE AND POTASSIUM CHLORIDE 
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Fic. 4. ASTM INDEX CARDS FOR NACL AND KCL 


may be applied in these cases, but, since the crystals 
within the metal sample may be preferentially ori- 
ented, care must be taken in comparing intensities 
with those in the ASTM index. 
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Fic. 5. X-RAY DIFFRACTION PATTERN OF PUGET SOUND GLACIAL CLAY 


TABLE I 
INDEXING OF X-RAY DIFFRACTION PATTERN 
(Figure 3) 
d-SPACINGS FROM ASTM 
OBSERVED DATA Innex 
20 d NaCl KCl 
27.3 3.26 7 3.25 aieey: 
28.4 3.14 100 seat 3.14 
31.7 2.82 2.81 
40.5 2.23 48 owe 2.21 
45.4 1.99 40 1.99 ae 
50.2 1.82 11 BeAr 1.81 
53.8 1.70 1 1.70 eee 
56.4 1.63 10 1.63 oe 
58.6 1.37 6 1.57 
66.3 1.41 13 1.41 1.40 
73.6 1.29 5 1.29 1.28 
75.3 1.26 6 1.26 prs: 
83.9 1.15 4 1.15 a 
87.6 1.12 2 ae 1.11 


A result of this type of work is shown in Fig. 5, 
a pattern of a sample of Puget Sound glacial clay 
obtained by E.C.S. Rao, EES research assistant, in 
connection with his project (No. 111) on the funda- 
mental properties of this material. This example 
shows one of the limitations of the method when used 
alone, i.e., the fact that some patterns contain so 
many lines that another identification tool must be 
used to aid in the determination. In this case both 
polarizing and electron microscopy were utilized. 

It should be noted here that one of the chief ad- 
vantages of qualitative diffraction techniques is the 
ability to differentiate between isotropic phases pres- 
ent in a sample or to determine, in the case of ceramic 
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reaction, what crystalline phases are present when 
the reaction is complete. 
Quantitative Analysis 

Quantitative analysis by X-ray diffraction methods 
is a very precise procedure, providing the proper 
instrumentation and techniques are utilized. It is 
based on the discovery that the intensity of a given 
line is directly proportional to the concentration of 
the substance producing the line. This statement 
somewhat oversimplifies the exact nature of the prob- 
lem, but, for practical purposes, is true. 

A series of standard curves are prepared by meas- 
uring the ratio of the intensity of the line in question 
for varying concentrations of the substance with a 
line of an “internal standard” whose concentration 
is constant. The accuracy of this method involves 
(1) the use of pure materials for preparation of 
standards, (2) the use of uniform, finely divided 
powders, (3) intimate mixing of the components, 
(4) careful mounting of the specimens, (5) instru- 
ment stability and reproducibility of results from 
day to day, and (6) care in measuring intensities. 

Although a great deal of time is sometimes con- 
sumed in the preparation of the standard curve, the 
actual analysis of unknowns, after sample preparation, 
can be made in a few minutes. 


Uses of X-Ray Spectroscopy 
Qualitative Analysis 
Since von Hevesy’s discovery of hafnium by means 
of X-ray spectroscopy, this analytical tool has been 
found extremely valuable, owing to the small number 
of lines in the spectra and the simple relation of line 
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Fic. 3. X-RAY DIFFRACTION PATTERN OF SODIUM CHLORIDE AND POTASSIUM CHLORIDE 


based on the fact that each crystal has its own charac- 
teristic patterns of d spacings and, hence, of the posi- 
tion of the diffraction lines. In addition, the location 
and types of atoms within the various planes cause 
the intensity of the diffracted lines to vary. It is then 
apparent that each crystalline material will have a 
set location of all its diffraction lines for a given wave 
length, and that each line will have a given relative 
intensity. The sample is usually finely ground (—200 
mesh) in order that the particles will have a high 
degree of random orientation. 

This randomness is necessary in order to compare 
patterns of the same material, since preferred orienta- 
tion would cause the relative intensity of certain lines 
to be lesser or greater than normal. Figure 3 shows 
a typical diffraction pattern of a sample containing 
two unknowns obtained in the Norelco Geiger- 
counter diffractometer in the X-Ray Diffraction Lab- 
oratory of the School of Mineral Engineering. 
Indexing of this pattern results in the data shown in 
Table I. A set of index cards has been published by 
the American Society for Testing Materials which 
gives the results of various diffractionists on more 
than 2,500 different crystalline materials. Each mate- 
rial is indexed according to the d value of the three 
most intense lines. Therefore, by selecting the strong- 
est lines of the pattern in Fig. 3 and checking the 
cards in the ASTM index, we find that the combined 
patterns from cards of sodium chloride and potas- 
sium chloride (Fig. +) match the lines in Fig. 3. 

All samples to be studied may not be in powder 
form, eg., metal specimens. The same procedure 
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may be applied in these cases, but, since the crystals 
within the metal sample may be preferentially ori- 
ented, care must be taken in comparing intensities 
with those in the ASTM index. 
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A result of this type of work is shown in Fig. 5, 
a pattern of a sample of Puget Sound glacial clay 
obtained by E.C.S. Rao, EES research assistant, in 
connection with his project (No. 111) on the funda- 
mental properties of this material. This example 
shows one of the limitations of the method when used 
alone, i.e., the fact that some patterns contain so 
many lines that another identification tool must be 
used to aid in the determination. In this case both 
polarizing and electron microscopy were utilized. 

It should be noted here that one of the chief ad- 
vantages of qualitative diffraction techniques is the 
ability to differentiate between isotropic phases pres- 
ent in a sample or to determine, in the case of ceramic 
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reaction, what crystalline phases are present when 
the reaction is complete. 


Quantitative Analysis 

Quantitative analysis by X-ray diffraction methods 
is a very precise procedure, providing the proper 
instrumentation and techniques are utilized. It is 
based on the discovery that the intensity of a given 
line is directly proportional to the concentration of 
the substance producing the line. This statement 
somewhat oversimplifies the exact nature of the prob- 
lem, but, for practical purposes, is true. 

A series of standard curves are prepared by meas- 
uring the ratio of the intensity of the line in question 
for varying concentrations of the substance with a 
line of an “internal standard” whose concentration 
is constant. The accuracy of this method involves 
(1) the use of pure materials for preparation of 
standards, (2) the use of uniform, finely divided 
powders, (3) intimate mixing of the components, 
(4) careful mounting of the specimens, (5) instru- 
ment stability and reproducibility of results from 
day to day, and (6) care in measuring intensities. 

Although a great deal of time is sometimes con- 
sumed in the preparation of the standard curve, the 
actual analysis of unknowns, after sample preparation, 
can be made in a few minutes. 


Uses of X-Ray Spectroscopy 
Qualitative Analysis 
Since von Hevesy’s discovery of hafnium by means 
of X-ray spectroscopy, this analytical tool has been 
found extremely valuable, owing to the small number 
of lines in the spectra and the simple relation of line 
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Fic. 6. X-RAY SPECTROGRAPH FOR QUALITATIVE ANALYSIS OF ORE 
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Fic. 7. X-RAY SPECTROGRAPH FOR QUANTITATIVE 
ANALYSIS OF TUNGSTEN IN TOOL STEEL 


intensity to concentration. However, the necessity of 
using the sample as the target of an X-ray tube and 
the time involved in photographically recording the 
spectrum on film have precluded any great use of this 
method. 

With the advent of high-intensity X-ray tubes of 
the Coolidge type and efficient Geiger-counter spec- 
trometers for accurately detecting and measuring the 
intensity of diffraction lines, the use of the fluorescent 
spectrum has become an accepted technique. The 
qualitative analysis of materials merely requires the 
operator to calculate the wave lengths of the spec- 
trum obtained from the sample and then determine 
what materials emit X rays of such wave lengths. 
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This procedure is limited at present, however, to 
analyzing for elements of atomic number 22 (tita- 
nium) and higher, because elements of lower atomic 
number emit X rays of such long wave lengths that 
they are absorbed in the air-path between the sample 
and the Geiger tube. 

While routine analyses in metals, ores and chemi- 
cals prove valuable, the greatest advantage is ob- 
tained in analyzing for those substances which are 
difficult to identify by other means. For instance, 
Fig. 6 shows an X-ray spectrograph of an ore sent 
to the laboratory for analysis which was found to 
contain iron, manganese, columbium, neodymium, 
yttrium, cerium, lanthanum, titanium, strontium, 
thorium, tantalum, and nickel. The lines of molyb- 
denum and copper, which form a constant back- 
ground, are due to the presence of these elements in 
the target of the X-ray tube. In these days of great 
interest in the rare earths and some of the rarer 
metals, such as tantalum and columbium, this method 
has proved invaluable to the prospector. 
Quantitative Analysis 

The nearly linear relationship between concentra- 
tion and intensity of a given spectral line is affected 
only by the absorption effects of other materials in 
the sample. For this reason a series of calibration 
curves may be necessary for certain analyses if the 
composition of the material varies. If, however, the 
composition is fairly constant, the use of a single 
calibration curve will suffice. 

The determination of major constituents of a 
sample (from a few to 100 per cent) gives X-ray 
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Fic. 8. CALIBRATION CURVE FOR QUANTITATIVE 
DETERMINATION OF MANGANESE IN CAST STEEL 


spectroscopy one of its greatest advantages over 
optical spectroscopy, since the former can be used to 
obtain results accurate within a fraction of a percent 
in a matter of 30 to 60 seconds. Figure 7 shows a 
calibration curve obtained in the determination of 
tungsten in tool steels, and Fig. 8 gives a similar 
curve for the analysis of manganese in cast steel. The 
latter was prepared by R. T. Torgerson, EES re- 
search assistant, for his work on Project No. 139. 

The analysis of elements of lower concentration is 
also possible, although requiring more careful tech- 
nique and more time. The method of analysis is 
identical, and the range of this instrumentation avail- 
able in the School’s laboratory makes determination 
possible down to about 10 parts per million, depend- 
ing largely, of course, upon the element under con- 
sideration and the matrix of the sample. 

Figure 9 presents a calibration chart showing the 
determination of Fe,O, in glass sands. A glass was 
melted from a batch of which 40 per cent by weight 
was silica. The molten glass was formed in an ap- 
propriate mold and, after cooling, was used as the 
sample. Although the standard method of reporting 
iron content in sand is in percentage of Fe,O,, the 
iron content of the glass is of the order of 1 per cent 
per 100,000. A comparison of the intensity of two 
lines is shown in Fig. 10. The standard sample is one 
made up in the laboratory of pure materials with a 
known amount of Fe,O, added, and the sand sample 
is one whose iron content was determined by wet 
analysis. 

Figure 11 shows a curve for the determination of 
nickel in raw ores. This curve was obtained by com- 
paring the intensity with the average of a number 
of wet analyses on the various samples. The deviation 
from the average here was about 0.1 per cent, as 
compared to a 2 per cent deviation from the average 
by other methods of analysis. 

(Continued on page 28) 
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THE IMPORTANCE OF CROSSWIND FORCES 
ATTAINABLE FROM MOVING SYSTEMS 


S. EASTMAN 


Professor of Aeronautical Engineering 


Any moving wing system, 
including the screw propel- 
ler, can create crosswind 
forces that offer important 
advantages over those ob- 
tained simply by inclining 
the plane of wing motion. 
An adaptation of blade 
element theory will illus- 
trate these forces and show 
the importance of making 
full use of them. 

Consider the action of a 
pair of wing segments that have equal and opposite 
crosswind velocity, as shown in Fig. 1. The pair may 
represent a segment of a single wing, first at an 
instant during its downward motion, and later at an 
instant when it attains an equal upward motion; 
or it may represent corresponding segments of two 
identical wings of a given moving wing system, one 
of which is descending at the instant that the other 
is ascending. Thus, a study of the forces acting on 
the pair can be used to illustrate the average force 
acting on a segment of a propeller blade at two 
instants, one a half-revolution later than the other; 
or it can be used to illustrate instantaneous forces 
acting simultaneously on corresponding segments of 
a contrarotating propeller; or it can be applied to 
very different types of moving wing systems in- 
cluding bird flight. 

Local disturbance of the fluid approaching the 
segments will be neglected for this illustration. 
Therefore, the relative velocity, », has the same 
magnitude for both segments, the advance velocity, 
Ua, is the same for both, and the crosswind velocity, 
Ye, is equal but opposite. The mean lift and mean 
thrust of the pair can be expressed in terms of the 
lift and drag of the individual segments: 


L=— (257+), (1) 


Ub 2 Ub 


Va La—La Va DatDa 
Ub ( 2 )- Ub ( 2 ) (2) 
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RESOLUTION OF DRAG FORCES 
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LIFT FORCES 
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DESCENDING SEGMENT ASCENDING SEGMENT 
Fic. 1. FORCES ACTING ON A PAIR OF WING SEGMENTS 


Since the drag terms are relatively small, lift de- 
pends mainly on the mean segment lift, Lz, and 
thrust depends mainly on the difference in segment 
lift, Lr, as follows: 


(4) 
La=Li+Lr, (5) 
(6) 


Independent control of lift and of thrust is possible 
with two controls. One will increase or decrease the 
angle of attack of both segments equally. The other 
will increase or decrease their angle of attack dif- 
ferentially. The former will control L,, and the 
latter will control Lr. These important components 
of individual segment lift are illustrated by the 
dotted arrows in Fig. 1. 
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RESOLUTION OF 
SEGMENT LIFT 
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Fic. 2. LIMITS OF ATTAINABLE LIFT AND THRUST AT 
DIFFERENT RATIOS v,/Vq 


The diamond-shaped areas in Fig. 2 outline the: 


values of lift and thrust attainab e when maximum 
positive segment lift equals maximum negative seg- 
ment lift. Boundaries of these areas are obtained by 
assuming Dz=D,=0 in Eq (1) and (2), and allowing 
first Lg and then L, to vary from maximum positive 
to maximum negative lift while maintaining maxi- 
mum lift on the other segment. It is important to 
observe that, when the crosswind velocity equals the 
advance velocity, v./va=1 and equal values of lift 
and thrust are attainable, but that smaller values 
of v,/va favor lift, and larger values favor thrust. 
Since the active portion of the blades of a high- 
pitch screw propeller operates near v,/va=1, a con- 
trol mechanism could be designed that would en- 
able each pair of blade segments to create a cross- 
wind force as large as its thrust. Such a mechanism 
has been developed for the helicopter but, owing to 
the very high values of v./v, that are used, the 
relatively small crosswind forces are obscured by the 
effect of tilting the rotor. If a moving wing system 
employing very small values of v./va can be de- 
veloped, it is possible that wing surfaces, moving in 
a way that will produce the needed thrust, can 
create the much larger crosswind force that is re- 
quired for lift. Bird flight provides the most common 
illustration of such a moving wing system, but pos- 
sibly the conventional screw propeller can be 
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Fic. 3. LIFT AND DRAG FOR THE PAIR OBTAINED FROM 
SEGMENT LIFT AND DRAG CURVES 


altered so that it will gain most of the attainable 
advantages without adding to the losses that 
normally arise in the production of thrust. 

The influence of segment drag will be illustrated 
by assuming that the increase in drag with lift is 
proportional to the square of the change in lift. 
Letting M be the proportionality factor, using the 
notation shown in Fig. 3, and expressing all forces 
in terms of minimum drag, Do, gives the following 
expressions for segment drag: 


DatDa _ Lor | Lr¥ 


Di-Da 2M Loi Lr 


2 Do 
Substituting Eq (3), (4), (9), (10) in (1) and (2), 
aad Va Ve 2M Lox Lr 


The parabola used in Fig. 3 to approximate the 


actual variation of drag with lift is determined by 
the value of M. By definition, 
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Fic. 4. INFLUENCE OF CROSSWIND VELOCITY ON LIFT 
From Eq (14), using M=0.0004 and Lor /Lr=1. 


Dor/Do__ Dor Dr 


In terms of wing-section coefficients this ratio is 
7 Cpo, for which 0.05 x 0.008=0.0004 is suf- 
ficiently representative to be used for the illustra- 
tions that follow. 

The ratio of the lift of the pair to fixed wing lift is 
of interest because it shows the gain in lift due to 
crosswind velocity. When the wing motion is 
stopped, v.=0 in Fig. 1, and the fixed-wing condition 
is represented. In this special case, % is reduced to 
v, and fixed-wing lift can be written 


L=( Li. (13) 


When used with Eq (11), this gives 


Ub Lox Lr 
The solid line in Fig. 4 is a plot of this expression 
with L7=0, so that the influence of the drag term is 
eliminated. The dotted line is indicative of the 
maximum influence of the drag term, because the 
assumed values, Lo,/L,=1 and Lr/Dos=50, are 
larger than those to be expected. Apparently the 
influence of blade drag may cause an appreciable 
increase in lift when the crosswind velocity is high, 


M= 


12 


relative to the advance velocity. However, the increase 
is small throughout the important range (v,/va) <2, 
and will be neglected for the sake of simplicity. 

The ratio of the thrust of the pair to the fixed- 
wing value of minimum drag will be used to illustrate 
the variation in thrust with crosswind velocity. 
When the wing motion is stopped, segment drag 
becomes the fixed-wing value, 


D,=( and Dy= (+ Do. (15) 


The ratio of thrust to Do; is obtained from Eq (15) 
and (12): 


T ve Lr Lot Lr\¥ 
The solid line curves in Fig. 5 are plots of Eq (16) 
illustrating the special case- when Lo;=0.. The 
dotted curves show the effect of operating at a 
higher angle of attack, so that Lp,/Do=50.. The 
pair overcomes its own drag when 7=0, and it is 
important to note that this condition requires very 
little crosswind velocity. Even when Lr is only 
about 1/15 of maximum segment lift (see curve for 
L7r/Do=10), the pair overcomes its own drag when 
the crosswind velocity is little more than 1/10 of the 
advance velocity. Birds employ a relatively slow 
wing motion while in steady horizontal flight. 
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Fic. 6. INFLUENCE OF CROSSWIND VELOCITY ON DRAG POWER 
From Eq (20) with Lop =0 and Mr=0.0004. 
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Probably they instinctively choose the method of 
flight that requires minimum propulsive effort. At 
any rate, the following study of segment drag power 
shows that minimum power for horizontal flight 
occurs at remarkably small values of v,/va. 

The mean segment drag power of the pair will be 
written, using Eq (9), as 


(17) 


Power required to overcome fixed-wing minimum 
drag is 


3 


Therefore, the ratio of drag power of the pair to 
fixed-wing minimum drag power is 


ao 


It may be possible to introduce periodic changes in 
the drag curve illustrated in Fig. 3 so that (Da+D,.)/2 
=D,, and Lr causes no change in drag. Birds may 
be able to do this because they have some control 
over the camber of their wings. A rigid wing with 
a spanwise hinge could accomplish the purpose 
nearly as well if the periodic change is small. At 
any rate, it appears to be possible to reduce, if not 
eliminate, the increase in drag with Lr. To represent 
such a condition the value of M should be reduced in 
the last term of Eq (19) but not in the middle term. 
Therefore two different values of M will be used: 


The solid line curve in Fig. 6 shows the increase 
in minimum drag power as v,/V. is increased. It is a 
plot of Eq (20) when both the lift and the thrust 
terms are zero. Some lift may be produced, because 
Ly will be the segment lift corresponding to its 
minimum drag. But the thrust must be the small 
negative value indicated by the curve for Lr/D)=0 
in Fig. 5, unless Mr can be reduced to zero. The 
dotted curves have been drawn to show the influence 
of the thrust term when Mr=0.0004, but they 
represent the influence of the lift term equally well 
when M,=0.0004. The broken lines have been 
added with the aid of a plot such as Fig. 5, to illus- 
trate the variation of 7/Dos; corresponding to the 
dotted curves. Minimum points on the broken lines 
move to the left and approach the solid line as Mr is 
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reduced, but, neglecting the possibility of such re- 
duction, the curve for T/Do;=0 shows that power 
required to overcome segment drag reaches its 
minimum value when 7,/va is less than 0.25. Even 
when 7/Do; is increased from 0 to 40, so that thrust 
may have a magnitude comparable with that 
normally required of lift, minimum power occurs 
when 2,/v. is less than one. Yet the moving wing 
systems employed for aircraft operate mainly in the 
region beyond v,/v.=2. 

The significance of this information is shown by 
Fig. 7, whith compares the wing profile drag power 
for hypothetical moving wing systems in which all 
segments have the same crosswind velocity. Fixed- 
wing drag power, which increases with v,°, provides 
the basis for comparison. Curve A represents the 
power required to overcome the profile drag of this 
wing by means of a conventional propeller which 
contributes an additional power loss assumed to be 
15 per cent of the wing drag power. Curve B repre- 
sents a moving wing system that operates with 
v-/Va=2, regardless of forward speed. If its total 
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wing area is the same as that of the fixed wing, its 
minimum wing drag power would be increased in the 
ratio (v%/va)® according to Eq(20), but its lift would 
be increased in the ratio »/va according to Eq (14). 
Therefore, it is assumed that wing area will be re- 
duced in proportion to the ratio of lift, L;/L=v./%= 
1/2.236 when v,/va=2. Consequently, profile drag 
power computed from the ratio (1/2.236) (v/va)*, is 5 
times that for the fixed wing. Curve C represents a 
moving wing system that operates with v./v,=2 at 
speeds up to } 1 max, but, at speeds beyond that for 
which increased lift is needed, it operates with 
v-/Va=0.2. This enables the same reduction in wing 
area as was assumed for Curve B, and the wing drag 
power is the same for speeds up to $v, max. The 
transition from v,/va=2 to 0.2 is made by using 
gradually diminishing values of v,/va as the speed 
increases, so that the increase in lift according to 
Eq (14) is just enough to provide the same maximum 
lift as that at } v, max. At speeds beyond the transi- 
tion, v-/va can be reduced to 0.2, and the drag power, 
computed from the ratio (1/2.236) (v%/v,)’, is only 
0.475 times that for the fixed wing. Apparently, for- 
ward flight can be accomplished with substantially 
less power than is required for the fixed wing, if 
optimum use can be made of the thrust and cross- 
wind forces attainable from moving wing systems. 


Dotted curves have been added in Fig. 7 to 
represent the relatively high thrust condition, 
T/Doy=10. Assuming that propeller losses remain 
15 per cent of the total thrust power of the screw 
propeller, they will increase in the ratio of the 
thrust, (7+Do,;)/Doy=11. Therefore, propeller losses 
have been increased: from 0.15 to 1.65 times that for 
overcoming wing profile drag. The corresponding 
increase for Curve B is negligible for this general 
comparison, because with v,/va=2, Fig. 5 shows that 
T/Doz=10 when Lr/D,y is much less than 10, and 
when translated to Fig. 6 this condition is represented 
by a point adjacent to the solid line curve. When 
v,/Va is variable, as assumed for Curve C, minimum 
Pp/Pos estimated from the curve for T/Doy=10 in 
Fig. 6 is about 1.6, and it occurs when v,/va is ap- 
proximately 0.45. In this case, therefore, wing drag 
power is (1/2.236) x 1.6, or 0.73 times the fixed-wing 
drag power, and it may be possible to reduce it con- 
siderably below this figure if Mr can be reduced as 
suggested previously. Apparently the potential ad- 
vantage of the moving wing system with variable 
v,/Va, When compared with the screw-propelled 
fixed-wing, is greater under conditions of high 
thrust than it is in level flight. This advantage is 
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important for steep climbs and for level flight at 
speeds less than the stalling speed of the fixed wing. 
Birds provide striking demonstrations of this fact. 

The foregoing analysis illustrates the action of a 
single pair of wing segments or a group of segments 
operating under identical conditions. By estimating 
the single pair that best represents the composite 
action of all segments, the illustration can be used 
to judge actual moving wing systems. A true 
quantitative study, however, requires the summing 
up of the action of many pairs. The wing action of a 
bird, for example, can be studied as completely as 
desired by considering bands of pairs located at 
different distances from the wing root, each of which 
is composed of enough individual pairs to represent 
the changes in crosswind velocity and wing configur- 
ation that take place during a complete cycle of 
wing motion. Similarly, the crosswind forces attain- 
able from a screw-wing (a term used to describe a 
wing system of the contrarotating propeller type) 
can be studied by considering several ring-shaped 
bands of pairs located at different radial distances 
from the axis. In each band the crosswind force 
produced. by an individual pair located on the 
horizontal diameter will act in the vertical direction, 
while that of a pair located on the vertical diameter 
will act in the horizontal direction, and all others 
will have components in both directions. By intro- 
ducing a periodic variation in the blade angle similar 
to the cyclic pitch control for helicopters, it is 
possible to control the crosswind forces on all pairs 
in such a way that the total crosswind force will act 
in any desired direction. Moreover, it appears to be 
possible to superimpose such additional periodic 
changes in blade angle that the thrust of pairs 
located at opposite ends of the vertical diameter, 
and of the horizontal diameter, can be controlled 
differentially, thereby providing positive control of 
both pitching moment and yawing moment. It has 
been pointed out that maximum advantage can be 
gained from crosswind forces when the value of 
v/a is very low. Therefore, the wind chords of a 
screw-wing should be set nearly parallel to the axis 
of rotation, like a screw propeller designed for a very 
high advance ratio. This permits the use of large 
highly tapered wings to satisfy strength and stiffness 
requirements. 

It is important to study the effect of inclining the 
plane of motion of the pair illustrated in Fig. 1, 
because an increase in maximum lift is possible that 
will benefit flight at low forward speed. Also, 
knowledge of this effect is needed for consideration 
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of the stability of moving wing systems. For the 
present illustrative purposes, however, it is sufficient 
to note that the drag power of the pair increases with 
the mean value of v,’, and consequently its minimum 
value occurs under the conditions assumed. Small 
angles of inclination result in some gain in lift with 
relatively small increase in drag power, but this 
simply means that the illustration presented here 
conveys a slightly conservative idea of the optimum 
characteristics attainable from moving wing systems. 

Between the instant when a wing segment acts as 
the descending, and that when it acts as the ascend- 
ing segment of a pair, it must undergo a change in 
velocity, or angle of attack, or both, in order to ac- 
complish the dual purpose of producing thrust and 
crosswind force. It is important to determine the 
influence of the resulting unsteady flow condition on 
the action of the pair. T. C. Lin, while employed as 
a research engineer at the University of Washington, 
completed a theoretical study! with which this in- 
fluence can be evaluated in the case of a wing sub- 
jected to harmonic crosswind motion. His work 
shows that the mean lift and thrust are not much 
affected if the wing chord is less than 1/15 of the 
wave length of the path, but the instantaneous values 
of lift are influenced in a manner that may prove to 
be very important. Apparently an additional peri- 
odic force is created by the unsteady flow. Under 
certain conditions this force directly opposes the 


inertia force that arises from the wing motion. 
Perhaps natural flight would be impossible if it did 
not take advantage of this method of relieving the 
seemingly intolerable inertia forces. At any rate, 
delicate insect wings withstand an acceleration of 
more than 2000 g if the frequency is 200 cycles per 
second with one-inch stroke, and the relatively 
heavy wings of a bird withstand an acceleration of 
more than 240g if the frequency is 20 cycles per 
second with one-foot stroke. It seems probable, 
therefore, that wings with harmonic crosswind mo- 
tion can be designed so that most of the inertia force 
is balanced aerodynamically all along the wing. If 
mechanical difficulties can be satisfactorily solved, 
this type of motion will provide maximum improve- 
ment in forward flight characteristics, because it 
enables the crosswind force of all pairs to be used 
most effectively for the creation of the needed lift. 
The cycloidal type of motion, recommended for 
many years by the late Professor Emeritus F. K. 
Kirsten,’ may offer the best solution to the mechan- 
ical difficulties, but all methods of producing periodic 
crosswind motion need to be considered. The author 
discussed several alternative moving wing systems in 
a memorandum’ issued three years ago regarding the 
practicability of the cyclogiro type of moving wing 
system tested in the F. K. Kirsten Aeronautical 
Laboratory. 
(Continued on page 28) 


M. I. T.’s Dr. Sherwood Delivers Lectures to Chemical Engineers 


D* THOMAS K. SHERWOOD, professor of chemical engineering and former dean 
of engineering at the Massachusetts Institute of Technology, was a visiting lec- 
turer in the Department of Chemical Engineering during the first two weeks of 
November. Dr. Sherwood, one of the country’s outstanding chemical engineers, has 
been at M.I.T. since 1930. During this period he has authored or co-authored two 
reference books in chemical engineering and over sixty articles in various technical 
publications. He has served as consultant to numerous industries and to the U. S. 
Government and has played an active role in many professional societies such as the 
A.L.Ch.E. and the A.S.E.E. He is especially well known for his work in engineering 
education and his studies in the field of mass transfer. His many honors include the 
William H. Walker Award in chemical engineering and the U. S. Medal for Merit. 

During his campus visit, Dr. Sherwood delivered four departmental lectures and 
conferred with the chemical engineering faculty and graduate students on research 
topics. “The Mathematical Treatment and Mistreatment of Data” discussed various 
methods used to correlate data in chemical engineering problems, especially prob- 
lems where large numbers of varables were suggested; ey certain potential dangers 
for their misue were pointed out. 

The remaining three lectures, “The Relation Between Mass Transfer and Fluid 
Friction in Turbulent Flow,” “Basic Concepts in Mass Transfer,” and “Barriers to 
Mass Transfer,” formed an integrated series in Dr. Sherwood’s field of specialization. 
Recent concepts of fundamental theory were reviewed for mass transfer between 
phases in turbulent motion with emphasis on the close relationship between this type 
of transfer and heat and momentum transfer under similar conditions. 

In addition to lecturing before the Chemical Engineering Department, Dr. Sher- 
wood presented addresses before the local sections of the A.I.Ch.E. and the A.C.S. 
and before the Vancouver Chemical Engineering division of the C.I.C. 


Dr. Thomas K. Sherwood 
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THERMODYNAMIC PROPERTIES OF ANILINE * 


JOHN H. LIENHARD, Jr. 
Engineering Experiment Station 
Research Assistant 


JosePpH L. McCARTHY 
Professor of Chemical Engineering 


WILLIAM B. NorDQUIST 
Assistant Professor of Mechanical Engineering 


Although aniline has long 
been known as a chemical 
especially useful as a base 
for dyes, only during recent 
years has it been investi- 
gated as a fuel.!:?:3 To facili- 
tate further study of aniline 
in this direction, available 
thermodynamic data have 
been collected, and, by 
methods of extrapolation 
and interpolation, the prop- 
erties of pressure, volume, 
temperature, enthalpy, and entropy have been es- 
timated in the gas and liquid phases of aniline be- 
tween 20° F and 2000° F and at pressures up to 
10,000 psia. 


J. H. Lienhard, Jr. 


Critical Constants and Phase Change Points 


The critical pressure and temperature of aniline 
were experimentally determined by Guye and Mallet 
as 769 psia and 797.8° F. These values have been 
quoted, and presumably accepted, by the editors of 
the International Critical Tables,’ by Kobe and Lynn,*® 
and by R. R. Dreisbach.”? In the present work the 
critical temperature and pressure were estimated by 
Meissner’s method’ wh'ch indicated values® with’n 
114 per cent and 3 per cent, respectively, of those 
determined by Guye and Mallet. 

The triple point among gas, liquid, and solid 
phases was located at 19.2° F by A. G. Chowdri and 
F. C. Auluck of the University of New Delhi.’° The 
normal ‘boiling point reported by a number of in- 
vestigators is as follows: 363.90-94° F,"! 363.78-90° 
F,” 363.9° 363.07° 363.74° F, and 363.92° 
F.° The value 363.9° F is accepted for this work. 
The normal freezing point has also been determined 
several times as follows: 20.7° F,” 20.8° F,” and 
21.0° F."6 Parks, Huffman, and Barmore’s value of 


*Condensed from the thesis of John H. Lienhard, Jr., sub- 
mitted for the degree of M.S. in Mechanical Engineering. 
The project, E.E.S. No. 150, was supervised by Dr. McCarthy, 
under whose direction the thesis was prepared, and Professor 
Nordquist. 
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20.7° F is the most widely 
accepted of these. 

No direct measurements 
seem to have been made of 
the critical volume. Since 
the rate of change of volume 
with temperature becomes 
immense near the critical 
point, this constant is best 
determined not by experi- 
mental methods but by the 
use of the Cailletet-Mathias 
law of rectilinear diameters. 
Applying this method, Guye and Mallet‘ arrived at 
the value 0.0471 ft/lb. Dreisbach’s equation for the 
Cailletet-Mathias rectilinear diameter’ yields the 
value, 0.0472 ft?/lb. The writers, using available 
experimental data, have written the following equa- 
tion for this line: 


J. L. McCarthy 


pit py=65.90—0.0303T, (1) 


where p; and p, represent the density in liquid and 
gaseous phase, respectively. Solution of Eq (1) at 
the critical temperature gives the critical volume as 
0.0473 ft?/lb. Meissner’s method yielded the value 
0.0502 ft?/lb.2 Dreisbach’s value of 0.0472 ft*/Ib is 
accepted for the present work. 


Pressure, Volume, and Temperature Relationships 


Saturated Vapor and Saturated Liquid 

The first step in the evaluation of pressure, vol- 
ume, and temperature data was that of determining 
the properties of the saturated-vapor and liquid by 
initially determining the vapor pressure-temperature 
relationship. From data of a number of experi- 
logarithm vapor-pressure versus 
reciprocal-temperature graph was made on a large 
scale.2 Two empirical expressions by Dreisbach’ 
were plotted along with these data and found to give 
results within 10 per cent of experimental values. 

These data are represented here with the following 
three linear equations 
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from 540° F to the critical point: 


P= 5.9088 — 3813/7, (2) 
from 180° F to 540° F: 

logio P=6.5141 —4409.2/7, (3) 
from 19° F to 180° F: 

logiy P=8.1660— 5469/7, (4) 


where P=psia, and 7=°R. In general these equa- 
tions represent the measured pressures within about 
2 per cent except between about 150 to 200° F. In 
this range Eq (3) and (4) yield values which may be 
as much as 8 per cent high. The rather sharp curve 
of the vapor-pressure line in this range tends to make 
the straight line approximations less accurate. 

Saturated liquid and vapor data were now com- 
pleted with the addition of volumes over the temper- 
ature range 19° F to 798° F. Saturated liquid density 
up to 200° F was determined from an equation for 
liquid density as a function of temperature, given in 
the International Critical Tables,> and based upon 
Tyrer’s data." Use of this equation was permissible, 
even though the expression was meant to apply at 
atmospheric pressure, since the vapor pressures are 
subatmospheric in the temperature range of interest. 
Density at these lower pressures will not vary sig- 
nificantly from values at one atmosphere. The 
reciprocal of these densities yielded the desired 
volume data. 

Saturated vapor volumes from 19°F to 200°F 
were estimated by the ideal gas law, inasmuch as no 
significant deviation from ideal behavior should exist 
here where the vapor pressure at 200° F is about 
0.67 psia. Gurevich and Sigalovskaya have given 
experimental values for aniline vapor densities at 
32° F and 122° F* which lie roughly within 1 per 
cent of ideality. It appears likely that the discrep- 
ancy here represents experimental error and that the 
gas is virtually ideal. 

From 200° F to the critical point, the saturated 
liquid and vapor volumes were calculated together, 
with the use of Eq (1). The vapor density was first 
estimated by use of generalized compressibility 
charts,”!'*.*8 and in the present work Obert’s chart™ 
was used. Liquid density was estimated® by sub- 
stitution of the vapor density in Eq (1). 

The data for saturated aniline are thought to be 
correct within 2 per cent from the boiling point to 
750° F, and within 5 per cent from here to the 
critical point. Below 364° F, 1 per cent accuracy 
can be claimed, based upon the following spot check 
at the normal boiling point by use of the Clausius- 
Clapeyron equation: 


dP/dT =AH,/TAV. (5) 
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The rate of change of vapor 
pressure with boiling point, 
dP/dT, has been experiment- 
ally determined by a num- 
ber of 
Timmermans’ value" of 19.6 
mm Hg/°C being repre- 
sentative; normal latent 
heat of vaporization (AH,) 
has been reported by Loug- 
uinine,* Awbery and Grif- 
fiths,?” and Marshall. Miss 
Marshall’s value of 113 
cal/gm is the most widely accepted; the volume 
change on vaporization, AV, calculated from these 
numbers was 6.30 ft/lb, a result in good agreement 
with the volume change of 6.28 ft*/Ib computed by 
use of the generalized compressibility factors. 

To aid in the interpolation of saturated vapor 
data below 350° F, the following two equations were 
developed® by combining the ideal gas law with 
Eq (3) and (4): 


W. B. Nordquist 


__1.2729 (10)8 


Below 180° F: p, (6) 


_ 2.8367 (10)’ 


180° F to 350° F: p, 


(7) 
A logarithm—P versus reciprocal-temperature plot® 
provided this additional value: 

305° F to 570° F: login py =4.108—4032.3/T. (8) 
Equation (6) is believed to give density values ac- 
curate within roughly 1 per cent at temperatures 
below 150° F. Equation (7) is accurate to about 1 
per cent above, and to about 3 per cent below 400° F. 

The pressure, volume, and temperature relations 
for saturated aniline are given in Table I. 


Superheated Gas 

No experimental data could be found for gaseous 
aniline. Thus pressure-volume-temperature relations 
were estimated entirely by use of generalized com- 
pressibility factors. Results are included in Table IT. 
Pressure-volume-temperature re'ations were also es- 
timated® by using the Beattie-Bridgman equation of 
state with constants evaluated by the generalized 
method of Maron and Turnbull,” to give the follow- 
ng statement: 


p= :08026[', _ 2.562 (10)? ] [ v+.334 (1 .231 )] 


vr? V 
36.75 [,  .2278 
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where: P=atmospheres, V=liters/gm-mole, and 
T=°K. This expression was checked at representa- 
tive points and found to agree with generalized com- 
pressibility values within about 1 per cent. Near the 
critical point, however, accuracy became exceedingly 
poor. On the 798° F isotherm, for example, a 1 per 
cent disagreement existed up to 600 psia, then rapidly 
increased to 90 per cent at the critical pressure. 


Subcooled Liquid 


Density of liquid aniline was estimated® by use of 
Watson’s expansion factors.*® At low pressures these 
values were checked at some points with the previ- 
ously mentioned liquid density expression’ and with 
data given by Gibson and Loeffler.** Certain other 
properties of liquid aniline have been reported: 


viscosity, sonic velocity,®*® and surface ten- 
sion.*2:40,41,42 


Enthalpy and Entropy of Vaporization 


The latent heat or enthalpy of vaporization was 
estimated® by Watson’s latent heat equation, 


AH,=AH,° (1—-T/T,)™, (10) 


the constants for which have been determined by 
Dreisbach’ as follows: 


AH, =194.1 (1— 7/699) (10a) 


The Clausius-Clapeyron equation was used to 
check results from this expression over the range of 
temperatures from 19°F to 798°F. Use of the 
Clausius-Clapeyron equation, however, required 
dP/dT values. Dreisbach has presented the following 
equation for dP/dT as a function of vapor pressure: 


dP /dT = .001374 (7.24179—logyP)? . (11) 


Values of dP/dT, obtained by differentiation of 
Eq (2), (3), and (4), ranged from about 15 per cent 
low to 20 per cent high when compared with the 
results from use of Dreisbach’s equation. This wide 
variation was expected, since Eq (2), (3), and (4) 
are only straight line approximations of a curve and 
thus will slope too much at the high temperatures 
and too little at the low temperatures. One value of 
dP/dT, available from an experiment by Timmer- 
mans! made at atmospheric pressure, agreed within 
11% per cent that obtained by use of Eq (11). Con- 
sidering these factors, the best estimate of dP/dT 
was taken to be that given by Eq (11), and these 
values, together with volumes appropriate for satur- 
ated liquid and vapor phases as shown in Table I, 
were used to estimate enthalpy of vaporization.°® 
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When the results obtained by use of Watson’s 
equation were compared? with those obtained by the 
Clausius-Clapeyron equation, the maximum varia- 
tion proved to be about 4 per cent up to 650° F. 
Equation (10a) gives a value of 203.5 Btu/Ib at one 
atmosphere, which agrees with the value obtained 
experimentally by Marshall.” 

Entropies of vaporization, AS,, were computed 
from the enthalpies of vaporization by the relation- 
ship, 


AS,=AH,/T. (12) 
Results are shown in Table I. 


Enthalpy 

The enthalpy of aniline has been arbitrarily fixed 
at zero for the standard reference state taken as the 
saturated liquid at the triple point, i.e., at a standard 
state temperature, To, of 19° F!° and a standard state 
pressure, Po, of 0.000592 psia or 0.001206 in. of Hg, 
from Eq (4). For any other state, the enthalpy can 
be calculated by use of the equation: 


T 
To Po 


where the first term, AH,o, is the enthalpy of vapori- 
zation at the standard state pressure and is found by 
Eq (10a) to be 273 Btu/lb; the first integral gives 
the enthalpy of aniline at a particular temperature, 
but at a sufficiently low pressure so that its behavior 
is substantially that of an ideal gas; and the second 
integral adds the enthalpy contribution arising from 
increase in pressure above the region of ideal gas 
behavior. 

To evaluate the first integral, values for the specific 
heat at low pressure, Cy, in the gas phase are needed. 
For the temperature range of present interest, no 
experimental data seem to have been published. The 
necessary specific heat values have therefore been 
estimated by two separate methods: by use of 
generalized bond frequencies,“ and by considering 
contributions of specific chemical groups to specific 
heat.** The results of the two methods agree with 
each other within about 3 per cent over practically 
all of the temperature range of interest. Results from 
the former are thought to be the more accurate and 
have been accepted for this work. The actual 
integration of the first integral was performed 
graphically upon a plot of C, versus temperature.’® 

The second integral could be evaluated by per- 
forming the indicated differentiation and integra- 
tions graphically; however, accuracy of such a 
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TABLE I 
SATURATED VAPOR AND LIQUID ANILINE 
| ‘ 
UM 3 P -° 
Vo_LuME (ft?/Ib) ENTHALPY (Btu/Ib) Entropy (Btu/Ib-°R) 
°F Psia Vi A Vo Vo Hi AH» Hg AS» Sg 
20 000592 01532 93,500 93,500 0 273 273 0 .569 569 
30 00101 01542 900 55,900 5 272 277 .007 555 562 
40 00170 01547 33,900 33,900 9 270 279 015 .540 555 
50 00275 01554 21,400 21,400 13 268 281 550 
60 00444 01561 13,500 13,500 17 266.5 283.5 .031 513 544 
70 00684 01569 8,940 21.5 264.5 286 .040 499 539 
80 0108 01576 5,780 780 25.5 263 288.5 534 
90 0165 01582 3,830 3,830 30 261 291 054 .475 529 
100 01589 2,580 2,580 35.5 258 293.5 .461 525 
120 0547 .01606 1,200 1,200 45 255 079 .440 519 
140 1125 .01620 615 615 55 251 306 095 418 513 
160 . 200 .01638 325 325 63.5 247.5 311 .110 .399 509 
180 .389 ; .01654 190 190 72.5 244 316.5 .125 .381 506 
200 01671 115.4 115.4 82 240 322 .364 505 
250 2.02 01715 40.2 40.2 111 228.5 339 .183 321 
300 5.17 0175 16.3 16.3 137 218 357 .221 .288 509 
350 11.8 0180 7.68 7.70 170 , 376 . 260 .255 515 
363.9 14.7 0181 6.28 6.30 177 203 380 .271 246 517 
24.2 0185 3.86 3.88 200 194 395 .295 .227 522 
450 46.2 0191 2.11 2.13 232 181 413 .329 .199 530 
500 82.0 0197 1.23 1.25 265 169 433 364 .176 540 
550 130 0203 .750 .770 295 153 .397 .152 549 
600 206 0217 470 492 325 137 462 429 .129 558 
650 300 0230 . 296 .319 356 117 473 463 .106 569 
700 420 0247 .187 .212 389.5 96.5 4 .498 .083 581 
750 590 0277 0913 .125 428 65.5 493.5 .538 054 592 | 
798.1 769 0472 0 .0472 487 0 487 .586 0 586 
TABLE II 
SUPERHEATED GASEOUS ANILINE 
PREssURE | TEMPERATURE °F 
=x 
In. Hg] Psia | © | 250 | 300 | 350 | 400 | 450 | 500 | 550 | 600 | 650 | 700 | 750 | 800 | 850 | 900 | 950 | 1000/ 1500 | 2000 
V | 166.4] 178.1] 190.1] 201.5] 213] 237; 248] 260] 272] 283.8] 295.5 | 307.5 | 318.6| 330.2 | 342.7} 460] 576 
1 | 0.49 H 339| 357| 396| 417] 439] 462] 486] 510| 535] 588| 615| 643] 671| 700] 1013| 1362 
.561| .585| .609] .633| .657| .680] .703| .725| .747| .769] .791| .812] .833] .854| .874|.1.059] 1.217 
35.5| 37.8| 40.1| 42.4] 44.8] 47.2] 49.4] 51.7] 54.1] 56.4] 58.8] 61.1] 63.8] 66.1] 68.5 92| 115.2 
5 2.45|{ H 357| 376| 417| 439] 462| 486] 510] 535| S61] 5 615| 643} 671| 700] 1013] 1362 
.551| .575| .599| .623 669 1] .713| .735| .757] .778| .799| .820] .840| 1.025 | 1.183 
(Vv 17.7| 18.9] 20.15] 21.1] 22.25] 23.45] 24.7] 25.86] 27.1] 28.1] 29.4] 30.6] 31.8] 33.0| 34.2 46| 57.6 
10 4.91 H 357| 376| 396| 417] 439] 462] 486] 510] 535] 5S61| 588] 613| 643] 671] 700| 1013] 1362 
511] .535| .560] .584| .608| .631| .654] .676| .698| .720| .742| .763| .784] .805| .825| 1.011] 1.168 
Vv 12.52 | 13.30} 14.1] 14.84] 15.63 | 16.40] 17.16] 17.92] 18.70] 19.5] 20.2] 21.1] 22.0] 22.8] 30.6] 38.4 
15 7.37 H 376| 396| 417] 439] 462| 486] 510| 535] 561| 588| 643| 700] 1013| 1362 
525| .550| .574| .598| .621] .645| .667] .689| .711| .733] .754] .775| .796] .816| 1.001 | 1.155 
(Vv 9.36| 9.91] 10.5 | 11.08] 11.67] 12.30] 12.80] 13.46 | 14.02] 14.62] 15.2] 15.82] 16.42| 17.03] 23.0] 28.9 
20 9.82 H 376| 396| 417] 462| 485] 510] 535] S61| S588] 615| 643| 671] 700] 1013| 1362 
518| .544| .568] .592| .615| .638| .661| .683| .705| .727| .748| .769| .790| .810| .995| 1.153 
7.90} 8.35] 8.81] 9.29] 9.84] 10.34] 10.71] 11.18] 11.62] 12.0] 12.64] 13.12 | 13.62] 18.35] 23.1 
25 | 12.28], H 396| 417| 439 2} 485] 510) 561| 587| 614] 643] 671| 1013] 1362 
538| .563| .587| .610| .633] .655| .677| .700| .722| .743| .764| .785| .805 1.148 
¥ 6.58| 7.00] 7.33] 7.72] 8.21] 8.55] 8.95] 9.34] 9.72] 10.2] 10.57| 10.93 | 11.43 | 15.32] 19.23 
29.9| 14.7 H 396} 416| 439] 462] 485| 510] S61} 587] 614] 642] 670] 699] 1013] 1362 
———'——_|| S$ 534| .558| .583| .606| .629| .651| .673| .696] .718| .739] .760] .781| .801| .986] 1.144 
Psia 
{Vv 4.80| 5.18] 5.36] 5.65| 5.98] 6.26| 6.59] 6.84] 7.13| 7.48] 7.74] 8.06] 8.40] 11.25] 14.13 
20 H 395| 415| 438] 461| 484 533] 5 587| 614] 642| 670| 699] 1013] 1362 
| 527| .552| .576| .600] .623| .645| .667| .689| .711| .732] .754| .775| .795| .980] 1.138 
( Vv 3.33| 3.52| 3.71] 3.96] 4.15] 4.36] 4.55] 4.74] .496| 5.15| 5.37] 5.59] 7.51] 9.42 
30 H 437| 460] 484] 509| 531] 560| 614] 642] 670] 699| 1012| 1361 
542| .566| .590| .613| .636| .658| .680| .702| .723| .744] .765| .785| .971] 1.129 
(Vv 2.48} 2.61] 2.75| 2.96] 3.10| 3.26| 3.39] 3.54] 3.70] 3.85] 4.05] 4.18] 5.62] 7.06 
40 H 413| 436] 458| 507] 530] 586| 613] 641| 669| 699] 1012] 1361 
534| .559| .5 7| .629] .652| .674| .696| .717| .738| .759| .779| .965| 1.123 
(Vv 2.07 | 2.18] 2.355| 2.45] 2.59] 2.70] 2.85] 2.95] 3.07] 3.20] 3.34] 4.49] 5.65 
50 { H 435| 529| 558] 585] 613| 641] 669] 1012] 1361 
553| .578| .601] .624| .646| .669| .691| .712| .733] .754| .774| .960] 1.118 
(Vv 1.702 | 1.793 | 1.930] 2.02] 2.14] 2.24] 2.34] 2.45] 2.54] 2.66] 2.77| 3.74] 4.71 
60 { H 456] 481] 504| 527| 557] 585| 612| 640| 669| 698] 1012] 1361 
s .573| .597| .619| .642|] .664| .687| .708| .729| .748| .770| .955| 1.114 
Vv 1.441 | 1.520} 1.648] 1.72] 1.833] 1.91] 1.995] 2.09] 2.17] 2.27] 2.37] 3.21] 4.03 
70 { H 433| 455 526| 556| 584] 61 698 | 1012] 1361 
.569| .593] .616| .638| .661] .683| .705| .726| .747| .767| .952] 1.111 
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TABLE II — Continued 


SUPERHEATED GASEOUS ANILINE 


PressurE| 5 TEMPERATURE °F 

PsIa A | 250 | 300 | 350 450 | 500 | 550 | 600 | 650 | 700 | 750 | 800 | 850 | 900 | 950 | 1000 | 1506 | 2000 

1.243 | 1.310] 1.410] 1.492 | 1.587 | 1.653 | 1.732} 1.82] 1.893] 2.26] 2.07] 2.80] 3.53 

80 H 432| 453| 479] SO1| 554] 584] 611] 640] 697] 1361 

5 .565| .590} .613| .635| .658] .680| .701| .722| .747| .764| .949) 1.129 

4 1.150 | 1.248 | 1.319 | 1.420] 1.462 | 1.534] 1.61] 1.675| 1.75] 1.83] 2.48] 3.14 

90 H 452| 477 523| 552| 583| 611| 639| 668} 697] 1011! 1361 

.563| .587| .610| .633| .655| .678| .699| .720| .744| .761| .947] 1.106 

1.023 | 1.115 | 1.174] 1.257 | 1.310 | 1.373 | 1.444] 1.500 | 1.572] 1.638] 2.23] 2.82 

100 H 451| 499| 522] 550} 583| 610] 639] 667] 1011 | 1361 

.560| .584| .608| .631| .653| .676| .697| .718| .739| .759| .945| 1.104 

Vv .836| .910| .969] 1.036| 1.079 | 1.148 | 1.197 | 1.244] 1.302] 1.360] 1.86] 1.76 

120 H 449| 473| 497| 521] 549| 582| 609] 638| 666] 696) 1011} 1360 

.554| .579| .602| .626| .649| .671| .693| .713| .735| .755]| .941] 1.099 

.766| .816| .874] .920| .961]| 1.018] 1.060 | 1.107] 1.161] 1.59] 1.562 

140 H 471| 495] 518| 547] 581| 608| 637] 666] 1010| 1360 

.574| .598| .622! .645| .668] .689| .710| .732] .752] 1.096 

.654| .698| .752| .796| .833| .896] .920| .962] 1.008] 1.39] 1.406 

160 H 468| 492| 515| 544| 579] 608| 636| 694] 1010} 1360 

.569| .594|] .618| .641| .664] .686| .707| .726] .748| .935| 1.093 

V .506| .610| .661| .700| .732] .772| .811| .849] .896| 1.23] 1.120 

180 H 462} 489] 512| 542| 607| 664] 693] 1009/ 1359 

.560| .590| .614} .638| .661| .682| .703| .725| .745| .932| 1.088 

Vv .543| .558| .620| .653| .690| .723| .760| .794] 1.106] .933 

200 H 487| 510| 541| 577| 605] 634| 693] 1009] 1358 

587} .611| .635| .658| .680| .701| .723| .743| .929) 1.083 

Vv .414| .447| .479] .508| .539] .569] .596| .628| .879| .799 

250 H 481 536| 574| 603} 661| 691] 1008] 1357 

s .579| .604| .629| .652| .673| .695| .717] .735| .924]| 1.076 

(Vv .325| .356| .387| .411| .440] .464] .515| .731| .698 

300 H 475| 499] 531| 572| 600] 630| 659} 689] 1007] 1357 

.569| .597| .622| .646| .668| .690| .712| .732] .916]| 1.072 

.290| .316] .341] .367| .389| .410] .434|] .622] .620 

350 H 493| 527| 567| 597| 657] 687! 1006] 1356 

s .591°| .617| .642] .664| .686| .708] .728| .913]| 1.070 

Vv .234| .262| .289] .310| .333| .352] .373| .542| .556 

400 H 488| 522|} 563| 594| 625] 655| 687] 1005} 1355 

s .584| .612| .637| .660| .682| .704| .724| .911| 1.067 

"4 .220| .245| .268| .288| .307| .326| .479| .S05 

450 H 516} 559] 590| 622| 652] 683] 1004] 1355 

s .607 | .633| .656| .679| .701| .722| .907| 1.065 

185] .209| .232] .250| .270| .288| .429|] .556 

500 H 508| 554] 587| 619| 650] 680] 1003] 1355 

.598| .628| .651| .674| .697| .717| .907]| 1.067 

.152| .1793| .203| .222| .240] .255] .389] 

550 H 501| 549| 583] 616] 647| 678] 1001] 1355 

.593| .624| .647| .671| .694] .715]| .905| 1.065 

1537| .180] .201| .215| .230] .354| .462 

600 H 543] 578| 644| 1001| 1354 

s .620| .644| .668| .691] .712| .904] 1.063 

1306| .158] .177] .193| .209] .326| .425 

650 H 532| 573] 608| 642] 673| 1000| 1353 

.640| .664] .688| .709] .901 | 1.061 

Vv 1094] .140| .159]| .1747] .1893| .300] .394 

700 H 524| 605| 639} 670] 998] 1353 

s .608| .635| .661| .685| .705] .899]| 1.059 

0472} .119] .139] .1538] .1674] .272| .359 

769 H 487| 598! 634| 667] 997] 1352 

.586| .629| .656| .681| .702| .897]| 1.057 

Vv .0378 | .1087 | .1293 | .1451| .1582] .261] .343 

800 H 466} 553] 595| 665| 997] 1351 

561| .626| .654| .679| .699] .895) 1.056 

0339 | .0771 | .1052 | .1212]| .1345] .228] .305 

900 H 452} 529] 584] 625| 659| 995] 1350 

.554| .607| .646] .673| .696| .893] 1.053 

0325 | .0554 | .0854| .1018| .1155| .205| .273 

1000 H 449| 571| 652] 993] 1349 

.552| .591| .632| .665| .691] .890| 1.051 

id .0285 | .0305 | .0332 | .0374 | .0423 | .0958 | .1344 

2000 H 433| 480] 522] 605| 975] 1339 

.539| .569| .597| .626] .651] .868]| 1.032 

Vv 0250 | .0261 | .0270 | .0281 | .0301 | .0500 | .0694 

4000 H 430| 474] Si1| 547| 585] 950) 1327 

s .529| .555| .575| .608| .632| 1.013 

Vv 0234 | .0240 | .0248 | .0257 | .0269 | .0387 | .0510 

6000 H 473} 545| S581} 942] 1318 

.519| .547| .572| .600| .626| .831 1.001 

lf v .0222 | .0231 | .0239 | .0250 | .0256 | .0334 | .0425 

8000 H 437| 514] 549| 584} 941] 1317 

518| .544| .572] .600| .624| .825| .995 

"4 0219 | .0223 | .0229 | .0236 | .0244| .0311 | .0377 

10000 H 441| 482] 517| 554] 589] 946] 1320 

514| .540| .568] .595| .618| .818]| .987 
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solution was questionable. Also, since the pressure- 
volume-temperature data had been determined by 
generalized compressibility factors, and since charts 
of pressure-correction versus reduced pressure for 
“general” gases exist in the literature, such a solution 
is unnecessary. The best chart of this kind found 
was one given by York and Weber.* In this chart 
(AH)M/T® was plotted against P, to give a family 
of constant 7, lines. Here, AH represents the re- 
duction of enthalpy due to change in pressure, 
M converts this change to a molar basis, and @ is a 
correction found necessary to make the data general- 
ize properly: @=(7./370)". Values for the exponent 
n were determined empirically on the basis of a large 
number of compounds. No values, however, were 
given below the critical point. Thus enthalpies down 
to the 800° F isotherm were established.’ Below the 
critical temperatures no ® exists, necessitating the 
following extrapolation: The 800°F enthalpy iso- 
therm was drawn in on a pressure-enthalpy plot.® 
Here was located the critical enthalpy. At the triple 
point pressure, the zero pressure enthalpies were 
marked in. Since the enthalpy did not vary more 
than 1 Btu/lb below the saturation pressure until a 
boiling point of about 400° F was reached, the iso- 
therms could be drawn in straight vertical lines 
below this temperature. The points where these lines 
run into the saturation line were also determined 
since, at the temperature of each, the vapor pressure 
was known. These points could now be connected 
to form the lower portion of the saturated vapor line. 
From here, with knowledge of the approximate shape 
that this line had to assume from other substances, 
and the critical point, the complete saturated vapor 
line could be interpolated. As the latent heats of 
vaporization had also been determined, the satur- 
ated line could next be drawn. Finally the remaining 
enthalpy isotherms which had been drawn vertically 
to about 10 psia, could be curved to meet the satur- 
ated vapor line at the proper pressure, a further 
interpolation. 


This treatment is thought, at the worst, to be ac- 
curate within 2 per cent of the values which would 
have resulted had accurate low-temperature ® values 
been available. 


Enthalpies in the liquid range were determined’ 
over a limited temperature range by use of another 
generalized method presented by Watson.*® Watson 
has generalized the difference between enthalpy at 
the given pressure and temperature, and that at the 
critical pressure and the given temperature. 
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Some measured values of C, have been reported 
for liquid aniline'’-*-®° at atmospheric pressure over 
the limited range of 0° F to 250° F. 


Entropy 
Entropy above the chosen standard state’ may be 
expressed by the equation, 


S= 


The first term was readily evaluated, since the 
standard state latent heat and temperature were 
known, i.e., 273/479=0.569 Btu/(lb) (°F). The 
second term was readily obtained by graphical in- 
tegration, Cy being previously evaluated over the 
complete range of temperatures. The third term, 
or the pressure correction, was evaluated by use of 
another generalized plot given by York and Weber.” 

Entropy of liquid aniline was also estimated for a 
few conditions’ by another of Watson’s expansion 
factor methods.*° 

Parks and Huffman* present entropy data at low 
temperatures and also a heat of formation for aniline. 


(14) 
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SEQUENTIAL WORK SAMPLING * 
BERL W. OWENS 
Assistant Professor of Mechanical Engineering 


Work sampling is a simple 
statistical technique used to 
evaluate how people use 
their time on their jobs. It 
is done in the following 
manner: First, the activities 
of a person, or group of 
persons engaged in similar 
work, are listed. For a sec- 
retary, the activities could 
be typing, filing, using the 
telephone, etc. An observer 
then looks at each of the 
subject workers at random times and enters a tally 
mark after the activity he sees being done at that 
instant. After many observations, the tally marks 
of each activity are counted. Time-utilization ratios 
are found by dividing the number of tally marks 
after each activity by the total number of marks. 

Utilization ratios furnish valuable information for 
management in measuring the effectiveness of man- 
agement or workers. For instance, a work sample 
taken of a group of production workers may show 
that delays are excessive. This fact may be obvious 
to someone familiar with the production operations, 
but it is very difficult for anyone to measure these 
delays from cursory observation. 

In the last three years the technique has gained 
widespread use, including that of sampling the ac- 
tivities of nurses and engineers, and promises to 
become one of the more valuable tools of the in- 
dustrial engineer. At the University, it is now 
included in the motion and time-study course. 


B. W. Owens 


Work Sampling 


The statistics involved in work sampling are those 
involved in general sampling theory. When we take 
many samples at random from a population of a 
known proportion, we expect the sample proportions 
to disagree with the population proportion in con- 
formity with binomial distribution. In work samp- 
ling we assume that the probability of observing a 
worker engaged in an activity is equal to the propor- 
tion of his time that the activity consumes. We 
specify that observations are to be made at random 

*From a paper presented by the author at the Regional 


Conference of the American Society for Quality Control, 
University of Washington, May 9, 1953. 
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times. Since utilization ratios are determined from a 
sample, the usual sampling errors are present, and 
hence we ask, “How reliable are these utilization 
ratios we have determined?” 

To answer this question from a statistical view- 
point, we apply the equation for the standard error 
of sampling from a proportion to the work-sampling 
situation: 


op=V/p (1—p)/N , (1) 


where 


o,=the standard error of a proportion, 
p=the probability of making an observation of 
an activity, 
N =the total number of observations made. 


As in any other random sampling situation, the 
size of the sample we take has a direct effect on our 
confidence that the sample is giving reliable informa- 
tion. In the case of work sampling we assume that 
the lot size is infinitely large. Considering the time 
utilization ratios determined by a work sample, we 
may say that we have a 90 per cent confidence that 
these ratios represent the true utilization ratios of 
the workers within some variation or tolerance that 
we can calculate. This plus or minus tolerance may 
be represented by the symbol B. The 90 per cent 
confidence level is equivalent to 1.645 standard 
errors of a normal distribution. We may determine 
a value for B as equal to 1.645 standard errors when 
using a sample size of N total observations and a 
proportion equal to the utilization ratio of an ac- 
tivity. The following symbols are used: 


n=the number of observations made of an ac- 
tivity with a subscript denoting a particular 
activity, b 

N =the total number of observations made for all 
activities, or the sample size, 

R=the utilization ratio for an activity=n/N, 

B=the ratio tolerance, plus or minus from R. 


Then, since 
1.645 
B=+/2.71R (1—R)/N. (2) 
To illustrate the use of the equations thus far 


presented, suppose that a work sample has been 
taken, that 280 observations have been made of a 
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certa‘n activity, and that a total of 2,000 observa- 
tions have been made of all activities. Then 
R=n/N =280/2,000 or 0.14. Applying Eq (2) for 
the tolerance, then B=+/2.71 (14) (.86)/2000 = 
0.000163 =0.0128. 

The utilization ratio may now be stated as 
R=0.14+0.013. This means that we have a 90 per 
cent confidence that the true ratio for the worker 
was somewhere between 0.127 and 0.153, but most 
likely, 0.14. If the value of N had been smaller, the 
value of B would have been larger. 

Whether or not we are concerned with this amount 
of variation depends upon what action may be taken 
by management as a result of having this information. 

The necessity of taking large samples to obtain 
acceptable tolerances is one of the greatest blocks to 
more extensive use of work sampling. To illustrate 
this problem more clearly, suppose we are sampling 
the work of 10 people engaged in similar activ.ty and 
we make 10 random observations of each during each 
day. We will collect 100 observations per day for 
the sample. In a period of a normal work week we 
will have 500 observations and in a month, approxi- 
mately 2,000 observations. A sample size of 2,000 
gives reasonable tolerances, but the ratios determined 
represent the activities of a month with little or no 
reliable information about daily or weekly variations 
in any one activity. Management must respect this 
possible error in acting upon work-sampling in- 
formation. With 40 subjects instead of 10, we could 
get the same sample size in a week and perhaps show 
significant differences in ratios from week to week, 
but, because few organizations have 40 workers in 
similar activities and the sample does not distinguish 
between individuals, this result also may be unsatis- 
factory. The application of work sampling may be 
greatly extended when fewer total observations are 
required for a given tolerance and when observations 
may be taken much more frequently than about 10 
per day per person. 

Sequential Work Sampling 

Sequential work sampling overcomes these disad- 
vantages... Under this plan, observations are made 
at regular intervals, let. us say, 5 minutes apart, 
instead of at random times. This sampling interval 
will yield 96 observations per person sampled per 
8-hour day. Now to ascertain the error and bias of 
a sequential work sample. 

Thorough examination of a worker’s activities will 
reveal two characteristics of importance to the 
statistical analysis of work sampling. First, the 
separate events of work that are compounded to 
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make up an activity begin at random times during a 
regular sampling interval if the interval is short 
enough, as illustrated by the situation of the tele- 
phone ringing on your desk. Although you may not 
receive telephone calls at random times during the 
day, any one telephone call will begin at a random 
time within, say, a 5-minute interval of time during 
the day. A telephone call will vary in length so that 
it will end at some random time within the same 
interval or one subsequent. In the light of this 
characteristic, it seems possible to sample at regular 
intervals and yet obtain a random, unbiased sample. 

The second important characteristic of worker ac- 
tivity lies in the way it is classified. Activities are 
distinguished by two general types: basic activities 
and interrupter activities. A simple illustration will 
show the difference between these two types. Con- 
sider a secretary doing normal office work including 
answering the telephone. She may type on a letter 
for several minutes at a time and stop only to 
answer the telephone. Activities such as typing may 
be called basic activities and those similar to answer- 
ing the telephone may be called interrupter activities. 
The amount of time used for an event of a basic ac- 
tivity is largely dependent upon the amount of 
interruption and, hence, is dependent upon the 
nature of the job. The amount of time used for an 
event of an interrupter activity is nearly independent 
of the over-all nature of the job and is closely related 
to the nature of the activity or the nature of an 
individual event of the activity. This second char- 
acteristic has statistical importance because it deter- 
mines the probability of making an observation of 
an event of an activity and thus determines the 
probability of observing the activity. It can be 
illustrated by again considering the activities of a 
secretary. Assume we are making observations every 
5 minutes and that one of our listed activities is for 
telephone calls. Assuming that somehow we have 
determined that the average telephone call takes 4 
minutes, then the probability that we would observe 
the secretary in the act of a telephone call is not the 
proportion of her total time that she spends on the 
telephone activity but, instead, is equal to the 
proportion of the sampling interval that the tele- 
phone call consumes. In other words, the actual 
probability of observing this activity is an effective 
probability dependent upon the average time for 
an event of the activity. In the example given, the 
effective probability is 4/5 minutes or 0.80. 

The estimate of the accuracy of a utilization ratio 
for an activity is changed when we can determine 
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WORK SAMPLING RECORD 


Job Sampled______Secretary Operator: N.A.N File No 
Observer: B.W.O. Dates: Page of. 
Remarks. Interval: 5 min 
Activity 30 40 70 80 90 100 R p’ 

| 20 | .208 60 
Filing ATH 6 | .063 | .500 
Dictation 1 | 010 | 010 
Misc. Prod. oH i 36 375 528 
Telephone 10 104 104 
Inquiries 7 073 145 
Personal 13 135 615 
Idle & A.D. ii | 3 .031 031 

TOTALS 96 


the effective probability of observing the activity. 
We can calculate the effect of knowing the effective 
probability by applying the equation for standard 
error when sampling from a proportion. The follow- 
ing symbols are used: 
Y =the number of events that are compounded to 
make up an activity, 
R=the utilization ratio n/N, 
p’ =the effective probability of making an observa- 
tion of an event, 
P,,=the probability of making u observations from 
Y events having an effective probability of p’. 


Then, 


Y! ‘\s (1 — 
(1—p’) 


The standard error of this probability distribution is 
(1—p’). 
The number of times that the activity will be ob- 


served is m, and an estimate of from Y events is 
n=Yp’. Substituting this in the standard error 


equation yields 
on=Vn (1-9’). 


Observing the activity u times will give R=n/N, so 
we may say that n=NR. Also, (or) =(on)/N. Sub- 
stituting these equalities in the above equation for 


on, We have 
Nor=VNR (1-?’) , 
or=VR (1—p’)/N. 


For the 90 per cent sampling confidence level, the 
value of the tolerance, B, equals the value of 1.645¢p. 
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Fic. 1. WorK-SAMPLING RECORD OF A SECRETARIAL JOB 


This may be substituted in the above equation for B: 


p- gy 


Determining the Accuracy of Utilization Ratios 
In random work sampling, Eq (3) applies by 
letting 


p'=R. (4) 
In sequential work sampling, the time-study method 
of determining p’ uses the following equation: 
p’=i/m, (5) 
where 
7=average time for an event of an activity, 
m=sampling interval. 
For activities of the interrupter type, 7 may be 
determined by conventional stop-watch time studies. 
However, values of p’ determined from i must be 
treated as estimates only, in that the time for one 
worker to perform an event may be different from 
that required by another, in other surroundings. 

A better method for determining 9’ is available 
when the time for an event frequently exceeds the 
sampling interval. This situation is revealed during 
sequential sampling, when several tally marks in a 
row are recorded for the same activity, usually in- 
dicating that an event has consumed more than one 
sampling interval and that there is a concentration 
of that particular activity during that part of the day. 

As an illustration, a sequential work sample of a 
secretary’s activity was taken with a sampling 
interval of 5 minutes. Figure 1 shows the record of 
that sample. Notice that this record shows concen- 
tration of the typing activity at definite times during 
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35,000 
INSTRUCTIONS FOR USING DIAGRAM 
TO DETERMINE RATIO TOLERANCE, B 
|. FOR EACH ACTIVITY, DETERMINE R FROM R= /N, 
2. IN SEQUENTIAL SAMPLING, DETERMINE 9’. 
32,500 N (IN RANDOM SAMPLING p' =R) 
3. ENTER THE LOWER SCALE WITH THE VALUE OF R 
AND PROCEED VERTICALLY TO A POINT (1) 
CORRESPONDING TO THE N OF THE SAMPLE. 
30,000 4. FROM POINT (1) PROCEED HORIZONTALLY TO 
THE RIGHT TO A POINT (2) DIRECTLY ABOVE 
THE VALUE OF p' ON THE LOWER SCALE. 
5. DETERMINE B BY EVALUATING THE POSITION 
OF POINT (2) AMONG THE TOLERANCE LINES. 
27,500 SAMPLE PROBLEM 
ASSUME: n=280, N=2000 OBSERVATIONS, 
p =0.40. 280/2000 = 0.14 
N 1) FIND R = 0.14 ON THE LOWER SCALE. 
2) FOLLOW DOTTED LINES ON DIAGRAM TO B =0.0106 
25,000 
| 
& 
22,500 
N 
20,000 
N 
17,500 
TN | 
\ 
15,000 
12,500 4 
10,000 \ 
! 
3 \ 
7,500 
\\ 
2,500 
0.30 0.40 0.50 0.60 0.70 0.80 0.90 Lo 


UTILIZATION RATIO, R AND EFFECTIVE PROBABILITY, p' 


Fic. 2. WorK-SAMPLING GRAPH 


the day. In such a case, the effect of this concentra- 
tion on the effective probability, p’, can be deter- 
mined by using the average number of observations 
made in a row for that activity during the taking of 
the sample. In Fig. 1 the activity of typing shows 
20 observations made in 8 groups. Let the symbol 
@ (a bar) equal the average number of tally marks 
in a row. Then for typing, equals 20 observations 
divided by 8 groups, or 2.50 observations per group. 
The equation that expresses this relationship is 
a=n/G, 


n=number of marks for the activity, 
G=number of groups of marks, one being 
the minimum sized group. 

From formal mathematical expressions of proba- 
bility we can derive the following expression for the 
average number of like things that will be drawn in 
a row from a proportion :* 

a=the average number of things drawn in a row 
in random sampling from a proportion, 
p=the proportion of like things in the population 
from which the sample is drawn. 
The drawing of things other than what we are in- 
terested in has no effect on d. This expression is fit- 
ting to the sequential work-sampling situation where 
p’, or the effective probability of observing an event, 
is the proportion in question. Conversely, if we 
know the value of G@, we can determine the value of 
the proportion. From the two expressions, 
a=1/(1-p’) and da=n/G, 


where 


an expression for p’ can be stated as follows: 


p’=1—(G/n). (6) 

The value of p’ for the typing activity in our 

example is 1—(8/20) or 0.60. The value of R for 

this activity is 0.208. The value of B, using p’ and 
substituting in Eq (3), is 


Neglecting the effective probability information and 
treating the sample as purely random yields 


B=.| 2.71 


The confidence in the resulting ratio, considering the 
sample as purely random, is poorer by 0.0682 —0.0485 
or 0.0197. This may be expressed as 0.0197/0.0485 
or 40 per cent poorer confidence. We now have two 
ways of determining p’ in sequential work sampling: 

p’=i/m (by time-study), (5) 

p’=1-—(G/n) (from sample). (6) 


*See the author's thesis for the degree of M.S. in Mechanical 
Engineering, University of Washington Press, Seattle, 1953. 
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Equation (5), which is estimated before the 
sample is taken, is less reliable than Eq (6) where 
a post-sample estimate is made. However, the latter, 
which involves group counting, is of little help for 
interrupter activities where the events are usually 
short and occur at random. Note that concentrations 
of such activities are not apparent on the work- 


_ sampling record (Fig. 1), as illustrated by the tele- 


phone activity. 

It is also important to note that the value of 
p’ can never actually be less than R, even though 
the equations say so. This also is illustrated by 
the telephone activity which shows 10 marks and 
9 groups. Then from p’=1—(G/n), p’=1—(9/10) 
or 0.10. In cases where p’ is estimated to be less 
than R, therefore, it is necessary to use Eq (4), 
p’=R, which was used in random work sampling. 

With the increased reliability of the utilization 
ratios as determined by a sample, it now becomes 
practical to apply control charts to worker activity 
in certain areas. To get a significant distinction be- 
tween the ratios as compared from one day to an- 
other, small tolerances must result for the ratios 
determined by one day of sampling. This demands a 
large number of observations and a high effective 
probability. These conditions can be met by using 
the sequential work-sampling plan with short samp- 
ling intervals or by the random sampling plan with 
a large group of similar workers. Sequential work 
sampling is particularly appropriate for sampling 
the work of a few people engaged in like jobs. 
Work-Sampling Diagram 

For any given sampling confidence, the equations 
for determining the’ ratio tolerances can be solved 
graphically. This graphical expression of the rela- 
tionship of B to R, p’, and N is shown by Fig. 2, a 
work-sampling diagram for a sampling confidence of 
90 per cent. Instructions for using the diagram are 
given in the upper right-hand corner. Tolerances for 
random samplings are determined at point (1), and 
for sequential sampling at point (2). Two prelimin- 
ary calculations are required for using the diagram: 

(1) Determine the utilization ratios for each ac- 

tivity from the equation, R=n/N, 

(2) Determine the values of effective probability 

for each activity by the most appropriate of 
Eq (4), (5), and (6): 
p'=R; p’=i/m;* p’=1—(G/n). 

With this diagram, mistakes are less likely to occur, 

and with a little instruction a clerk can do the job. 


*Equation (5), the time-study method, is seldom used. 


27 


| 


NORTHWEST ROAD CONFERENCE 
TO BE HELD ON CAMPUS 


HE SEVENTH NORTHWEST CONFERENCE ON ROAD 

BUuILDING will be held in More Hall, on February 
17, 18, and 19, 1954. Subjects to be discussed for- 
mally, or by panels, include Oregon and Washington 
highway problems, highway financing, urban express- 
ways, county road needs, aerial surveys, swelling 
sub-grade soils, roadside development, asphaltic road 
surfacing and performance, the engineer’s notebook 
(which will include descriptions of several unusual 
engineering situations connected with road building), 
and concrete and timber structures. Speakers especi- 
ally qualified in their fields from Oregon, Washing- 
ton, California and the East will be on the program. 
The highway directors of both Washington and Ore- 
gon, Mr. W. A. Bugge and Mr. R. H. Baldock, will 
address the group. Speakers from the Civil Engi- 
neering Department of the University will include 
Prof. M. I. Ekse, who will discuss methods of apply- 
ing the WASHO road-test data to roads and streets 
in Washington, and Prof. R. G. Hennes, who will 
speak on the work of the Washington Council for 
Highway Research, with special reference to the 
current project : the allocation of highway costs. 

Registration will take place in More Hall from 
8 :30 to 10:00 A.M., Wednesday, February 17. The 
registration fee of $2.00 will include a copy of the 
Proceedings of the Conference, when published. Road 
Building Conferences in the past have attracted a 
broad cross section of interest. Over 350 engineers, 
contractors, suppliers, and public officials attended 
the 1952 Conference at the University of Washington. 
The 1953 Conference was held at Corvallis, Oregon. 


IMPORTANCE OF CROSSWIND FORCES 
(Continued from page 16) 


At the present time, moving wing systems gain little ad- 
vantage from the attainable crosswind forces. If future 
developments can be made to yield the maximum advantage 
that is possible, aircraft of this type may have performance 
superior in all respects to that obtainable with screw-pro- 
pelled wings. In any event, important improvements should 
be expected from efforts to gain this advantage. 
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ASPECTS OF X-RAY DIFFRACTION 
(Continued from page 9) 

This account has merely given a cross section of the 
analytical work undertaken at the Laboratory during the 
past year. Studies of deformation and orientation due to 
metal working, solid solution determination, particle-size 
analysis, and a host of other uses not described here, attest 
the versatility of X-ray diffraction techniques. While not 
providing a solution to all analytical problems, they can be a 
powerful tool to materially assist the scientist as well as 
those in production. 


THERMODYNAMICS OF ANILINE 
(Continued from page 22) 
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CORRECTION 
In the October, 1953, issue of The Trend, Prof. F. C. 
Smith should have been cited as author of the technical 
article, “Simplified Solution of Castigliano’s Theorem,” 
Bull. Mech. Div., ASEE, Dec., 1951, which was listed under 
the name of Prof. S. Sergev in the report of the civil engi- 
neering publications. 
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